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Abstract. We show here the global, in time, regularity of the three dimen- 
sional viscous Camassa-Holm (Navier-Stokes-alpha) equations. We also pro- 
vide estimates, in terms of the physical parameters of the equations, for the 
Hausdorff and fractal dimensions of their global attractor. In analogy with 
the Kolmogorov theory of turbulence, we define a small spatial scale, 1^, as 
the scale at which the balance occurs in the mean rates of nonlinear transport 
of energy and viscous dissipation of energy. Furthermore, we show that the 
number of degrees of freedom in the long-time behavior of the solutions to 
these equations is bounded from above by (L/l^)^ , where L is a typical largo 
spatial scale (e.g., the size of the domain). This estimate suggests that the 
Landau— Lifshitz classical theory of turbulence is suitable for interpreting the 
solutions of the NS— a equations. Hence, one may consider these equations as 
a closure model for the Reynolds averaged Navier— Stokes equations (NSE). We 
study this approach, further, in other related papers. Finally, we discuss the 
relation of the NS— a model to the NSE by proving a convergence theorem, 
that as the length scale «i tends to zero a subsequence of solutions of the 
NS— a equations converges to a weak solution of the three dimensional NSE. 
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1. Introduction 

Proving global regularity for the 3D Navier-Stokes equations (NSE) is one of 
the most challenging outstanding problems in nonlinear analysis. The main difh- 
culty in establishing this result lies in controlling certain norms of vorticity. More 
specifically, the vorticity stretching term in the 3D vorticity equation forms the 
main obstacle to achieving this control. 

In this paper we consider a similar partial differential equation, the so-called 
viscous Camassa-Holm, or Navier-Stokes-alpha (NS— a) equations. The inviscid 
NS— a equations (Euler— a) were introduced in [22] as a natural mathematical gen- 
eralization of the integrable inviscid ID Camassa-Holm equation discovered in [2] 
through a variational formulation. Our studies in [4]-[6] indicated that there is a 
connection between the solutions of the NS— a and turbulence. Specifically, the 
explicit steady analytical solution of the NS— a equations were found to compare 
successfully with empirical and numerical experimental data for mean velocity and 
Reynolds stresses for turbulent flows in pipes and channels. These comparisons 
led us to identify the NS— a equations with the Reynolds averaged Navier-Stokes 
equations. These comparisons also led us to suggest the NS— a equations could be 
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used as a closure model for the mean effects of subgrid excitations. Numerical tests 
that tend to justify this intuition were reported in [7]. 

An alternative more "physical" derivation for the inviscid NS— a equations 
(Euler— a), was introduced in [23] and [24] (see also [5]). This alternative deriva- 
tion was based on substituting in Hamilton's principle the decomposition of the 
Lagrangian fluid-parcel trajectory into its mean and fluctuating components. This 
was followed by truncating a Taylor series approximation and averaging at con- 
stant Lagrangian coordinate, before taking variations. A variant of this approach 
was also elaborated considerably in [29] . See also [30] for the geometry and analysis 
of the Euler— a equations. For more information and a brief guide to the previous 
literature specifically about the NS— a model, see paper [18]. The latter paper also 
discusses connections to standard concepts and scaling laws in turbulence model- 
ing, including the relationship of the NS— a model to large eddy simulation (LES) 
models. Results interpreting the NS— a model as an extension of scale similarity 
LES models of turbulence are reported in [15]. 

Equations similar to the NS— a equation, but with different dissipative terms, 
were considered previously in the theory of second grade fluids [16] and were treated 
recently in the mathematical literature [8], [9]. Second grade fluid models are de- 
rived from continuum mechanical principles of objectivity and material frame indif- 
ference, after which thermodynamic principles such as the Clausius-Duhem relation 
and stability of stationary equilibrium states are imposed that restrict the allowed 
values of the parameters in these models. In contrast, as mentioned earlier, the 
NS— a equation is derived by applying asymptotic expansions, Lagrangian means, 
and an assumption of isotropy of fluctuations in Hamilton's principle for an ideal 
incompressible fluid. Their different derivations also provide the different inter- 
pretations of the parameter ai, namely, as a flow regime quantity for the NS— a 
equation, and as a fixed material property for the second grade fluid. 

The aim of this paper is to establish the global regularity of solutions of the 
NS— a, subject to periodic boundary conditions. We also provide estimates of the 
fractal and Hausdorff dimensions of their global attractors. In particular, we iden- 
tify the dimension of the attractor with the number of degrees of freedom governing 
the permanent regime of these equations and find a remarkable compatibility be- 
tween these estimates and the number of degrees of freedom in turbulence a la 
Landau and Lifshitz [27]. This leads us to regard the NS— a equations as a suit- 
able closure model for turbulence, thought of as an averaged theory rather than an 
individual realization, cf. [4], [5], [6], [23] and [24]. Finally, we relate the solutions 
of the viscous Camassa-Holm (NS— a) equations to those of the 3D NSE as the 
length scale ai tends to zero. Specifically, we prove that a subsequence of solutions 
to the NS— a model converges as ai ^ to a weak solution of the 3D NSE. 
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— {uqU — a^Au) — vA{aQU — a^Au) 



2. Functional setting and preliminaries 

We consider the following viscous version of the three dimensional Camassa 
Holm equations in the periodic box O = [0, L]^: 

di 

- ux {V X {alu- alAu)) + —Vp = f , (la) 

Po 

V ■ u = (lb) 

u{x,0)^u"'{x), (Ic) 

p TT 

where — — ^af^\u\'^ — a1{u- Au) is the modified pressure, while tt is the pressure, 

Po Po 

v > is the constant viscosity and po > is a constant density. The function / 
is a given body forcing > and ai > are scale parameters. Notice ao is 
dimensionless while ai has units of length. Also observe that at the limit ao = 
l,ai = we obtain the three dimensional Navier-Stokes equations with periodic 
boundary conditions. 

For simplicity we will assume the forcing term to be time independent, i.e. 
f{x,t) = fix). 

From (1) one can easily see, after integration by parts, that 



■ / {a^u — aiAu)dx = / fda 
Jn Jq 



d_ 

On the other hand, because of the spatial periodicity of the solution, we have 
Audx = 0. As a result, we have ^ udx = fdx; that is, the mean of the 
solution is invariant provided the mean of the forcing term is zero. In this paper 
we will consider forcing terms and initial values with spatial means that are zero; 
i.e., we will assume u'^dx = J fdx = and hence udx = 0. 
Next, let us introduce some notation and background. 

(i) Let X be a linear subspace of integrable functions defined on the domain f2, 
we denote 



X:={ipeX: [ ip{x)dx = 0} 
Jn 



(ii) We denote V = {ip : ip is a vector valued trigonometric polynomial defined on 
f2, such that V • = and ip{x)dx = 0}, and let H and V be the closures 
of V in L^(f2)^ and in H^{Q)^ respectively; observe that H-^, the orthogonal 
complement of H in L'^{nf is {Vp : p G H^(f^)} (cf. [11] or [31]). 

(iii) We denote : L^(f2)^ H the orthogonal projection, usually referred 
as Leray projector, and by ^ = —PaA the Stokes operator with domain 
D{A) = (i?^(rj))^ nF. Notice that in the case of periodic boimdary condition 
A = —A\jj^^^ is a self-adjoint positive operator with compact inverse. Hence 
the space H has an orthonormal basis {wj}^^ of eigenfunctions of A, i.e. 
Awj = XjWj, with < Ai < A2 . . . < Aj ^ oo; in fact these eigenvalues have 
the form \k\^^ with k G Z^\{0}. 

(iv) We denote (•, •) the L^-inner product and by | • | the corresponding L^-norm. 
By virtue of Poincare inequality one can show that there is a constant c > 
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such that 

c\Aw\ < \\w\\h^ < c~'^\Aw\ for every w G D[A) 

and that 

c\A^^^w\ < \\w\\hi < c~^|A^/^ti;| for every w eV . 

Moreover, one can show that V = D{A^/^), (cf. [11] , and [31]). We denote 
{{■,■)) = {A^^"^-, A^/^-) and |1 • || = \A^/^ ■ \ the inner product and norm on 
V, respectively. Notice that, based on the above, the inner product ((•,•))' 
restricted to V, is equivalent to the inner product 

[u,v] = al{u,v) + al{{u,v)) for u,vgV, (2) 

provided ai > 0. 

Hereafter c will denote a generic scale invariant positive constant which is 
independent of the physical parameters in the equation, 
(v) Following the notation for the Navier Stokes equations we denote B{u, v) = 
Pcr[(u • V)w], and we set B{v)u — B{u, v) for every u,v €V. Notice that 

{B{u, v),iu) = —{B{u, w), v) for every u, w, w e F (3) 

We also denote B{u,v) = —Pcr{u x (V x v)) for every u,v G V. Using the 
identity 

3 

(6 • V)a + ^ ajVbj = -b x {V x a) + V(a • b) , 

one can easily show that 

{B{u,v),w) = {B{u,v),w) — {B{w,v),u) 

= {B{v)u- B*{v)u,w), (4) 
for every u,v,w & V . As a result we have 

B{u,v) ~ {B{v) — B*{v))u for every u,v € V . (5) 

In the next Lemma, we show that the bilinear operator B can be extended 
continuously to a larger class of functions. 

Lemma 1. (i) The operator A can be extended continuously to be defined on 
V = with values in V C H''^ such that 

{Au,v)v' = {A^/'^u,A^/\) = [ {Vu:Vv)dx 

Jn 

for every u,v gV. 

(ii) Similarly, the operator can be extended continuously to be defined on D{A) 
with values in D{Ay such that 

(A^u^v) ]j(^j^y = {Au,Av) , for every u,v E D{A) . 

(iii) The operator B can be extended continuously from V x V with values in V , 
and in particular it satisfies 

\{B{u,v),w)v,\ < c\unu\\'/'\\v\\\\w\\ 
\{B{u,v),w)v>\ < c||u||||t;|||«;|^/^||w||^/^ , 
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for every u,v,w G V. Moreover, 

{B{u,v),'w)v' = —{B{w,v),u)v' , for every u,v,iu G V , 
and in particular, 

{B{u, v), u)v' = for every u,v €V . 
(iv) Furthermore, we have 



{B{u,v),w)D{Ay 



< c\u\\\v\\ Wwl 



,1/2 



\Aw 



1/2 



for every u G H, v €V and w G D{A), and by symmetry we have 



{B{u,v),w) 



< c\\u\ 



,1/2 



\Au 



1/2, 



for every u G D{A),v G V and w £ H . 
(v) Also, 

|V2|L,||l/2 



{B{u,v),w)D{Ay 



\v\\Aw\ + \v\\\u\\\\w\ 



1/2 



\Aw\'^^) 



for every u €:V,v G H,w & D{A) . 
(vi) In addition, 

{B{u,v),w)v\ <c(||w||l/2||Awf/2|i,|||u;|| + |Au||i,||w,|V2||y;||l/2j ^ 

for every u £ D{A),v G H,w € V . 

Proof The proof of (i) can be found in [11] or in [31]. The proof of (ii) is a straight 
forward extension of that of (i). 

To prove (iii), let us first consider the case when u,v,w G V. Then we have 

\{B{u,v),w)v'\ = I / u X (V X v) ■ wdx\ 
Jn 

< c\\u\\l3\\Vv\\l4w\\l» ■ 
Recall the following Sobolev inequalities in 

< c\Ml2^\Mh^ ' (6a) 

< c\\ip\\%^\\ip\\]^? , and (6b) 
llv'lUe < Mm , for every ^ G H\n) . (6c) 

Then by the above inequalities we have: 

\{B{u,v),w)v'\ < c|u|^''^||'u||^''^||t;||||w|| . 

Moreover, it is clear that for u,v,w G V 

{B{u,v),w)v' = -{B{w,v),u)v' ■ 

Since V is dense in V wc conclude the proof of (iii). 

Let us now prove (iv). Again we consider first the case where u,v,w G V 



\{B{u,v),w)n{Ay 



[m X (V X v)] ■ wdx\ 



< c||u||i^2||Vt;||L2||w;||i,c» 



Recall Agmon's inequality in R^: 



|</?||l~ < C||</9||^i^||</9||^2 



(7) 
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The above gives 

\{B{u,v),w)niAy\<c\u\ \\v\\ \\w\nAw\\'^^ . 
To prove ( v) we again take u,v,iu gV and we use (4) to find 



\{B{u,v),w)D{Ay\ < I / ({u ■ W)v) ■ wdx\ + \ / {{wV)u)-vdx\ 

Jn Jn 

< I / {{u-V)w-v)dx\ + \\v\\LA\^u\\LA\w\\Loo 

Jn 

< c||u||j;,3||Vw||i:,6|w| + c|t;| ||u|| ||t«||i:,oo • 

By (6b-c) and (7) inequalities wc finish our proof. 
The proof of (vi) is similar to (v). From (4) wc have 

\{B{u,v),w)v'\ < \ {{u ■ V)v)wdx\ + \ {{w ■ V)u) ■ vdx\ 



< I j {{u-y)w) ■vdx\+c\\'w\\L3\\'Vu\\Le\v\ 



By (6a) and (7) inequaUties we finish our proof. 



□ 



Wc apply Pa to (1) and use the above notation to obtain the equivalent system 
of equations 

^(q!oU + alAu) + vA{al + alA)u + 

+ B{u, alu + a\Au) = P^f , (8a) 
u{0) = u'" . (8b) 

Alternatively, if we denote 

V = a^u + a\Au (9) 

the system (8) can be written as 

dv 

— +yAv + B{v)u-B*{v)u = Paf, (10a) 

u(0) = u'" . (10b) 

Wc will assume that P^f = f, otherwise we add the gradient part of / to the 
modified pressure and rename Pa-f by /. 

Definition 2 (Regular Solution). Let f H, and let T > 0. A function u S 
C{[0,T);V)nL^{[0,T);D{A)) with ^ e L'^{[0,Ty, H) is said to be a regular solu- 
tion to (8) in the interval [0,T) if it satisfies 



(^("o" + "i-^u), w)d(a)' + i^{A{alu + a\Au),w)D(Ay 

+ {B{u, alu + ajAu), w)o{Ay =(/, w) , (H) 

for every w e -0(^4) and for almost every t € [0, T). Moreover, u(0) = w*" in V . 
Here, the equation (11) is understood in the following sense: 
For every to,t £ [0, T) we have 
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{aQu(t) + aiAu{t),w) + v / {aQu{s) + aiAu{s),w)ds+ 

J to 

'to 



■ / {B{u{s),alu{s) + alAu{s)),w)D{Ayds ^ / {f,w)ds . (12) 

J to ^ to 



3. Global existence and uniqueness 

In this section we prove global existence and uniqueness of regular solutions to 
equation (8). 

Theorem 3 (Global existence and uniqueness). Let f & H and u'" G V. Then 
for anyT > the equation (8) has a unique regular solution u on [0,T). Moreover, 
this solution satisfies: 

(i) u&L^,{{Q,T\-H\n)). 

(ii) There are constants R^, for fc = 0, 1, 2, 3, which depend only on v, ao, a\ and 
f, but not on u'", such that 



linisup (aol^^u]^ + al\A^uf) = 



for fc = 0, 1, 2, 3. In particular, we have 

1 . \ \A-y^f? kn^i ^ . r \f? \f? , 



that is: 

^2 G^iy^ . r 1 11 G^u^ 
i?o<^mm|-,,-^| = -p, 

w/iere G = — ^^tj «s f/ie Grashoff number, and 7"^ = min < ^ , \ . 
Fwrt/iermore, 

,, z-'+T --2 3 \ 1/2 
limsup- / (a2||«(s)f + a2|Au(s)nds<:.Aii?2<l^ (14) 

/or all t > 0. 

Proof. We use the Galerkin procedure to prove global existence and to establish 
the necessary a priori estimates. 

Let {wj}j^^ be an orthonormal basis of H consisting of cigcnfunctions of the 
operator A. Denote Hm = span{wi,... ,Wm} and let Pm be the L^-orthogonal 
projection from H onto Hm- The Galerkin procedure for the equation (8) is the 
ordinary differential system. 

^{alum + ctlAum) + vA{alura + a\Aum) + 

+ PmB{um,alum + afAum) = Pmf , (15a) 

Um{0) = PmU'^ ■ (15b) 

Since the nonlinear term is quadratic in Um, then by the classical theory of ordinary 

differential equations, the system (15) has a unique solution for a short interval of 
time {—Tm,Tm)- Our goal is to show that the solutions of (15) remains finite for 
all positive times which implies that Tm = oo. 
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-estimates 

We take the inner product of (15) with Um and use (4) to obtain 

\^i^o\'^rnf + allium 11^) + Z/(q;o ||^ + Q!?|AUm|^) = {Pmf,U„ 



Notice that 



|(Pm/,«m)| < <! j^-l/J/IIll™ 



and by Young's inequahty we have 

I 4-1 f\'2 

+ ^al\Aum\'^ 



\{Pmf,Um)\ < 



2||,, l|2 



M-l/2j|2 M-lf|2 

Denoting by Ki = min{ ^ , ^ — }, from the above inequahties we get: 

j^{al\um? + aWWraf) + z^(ag||w™f + al\Aum?) < Ki . (16) 
By Poincare's inequality wc obtain 

and then by Gronwall's inequahty we reach 

al\um{t)f+al\\um{t)f < e-''^^*(a^|„„(o)|2 + a2||^^(o)||2)^ 

(1-e- 



^1 -^.Al^^ 



That is 

al\um{t)\' + al\\um{t)f < h ■■= aly-\' + a?||«™f + ^ . (17) 

H"^ -estimates 

Integrating (16) over the interval {t, t -\- t) 

rt+r 

V I (a^||w„(s)||2+a?|Au^(s)|2)ds<rifi + (a^|w„(t)|2 + af||u„(t)f) 

< tKi + ki hir) . (18) 
Now, take the inner product of (15) with Aum to obtain 

^^("ollwmf + al\Aumf) + u{al\Aum\^ + al\A^/'^Umf) + 

+ {B{Um, alum + a\AUm)-, AUm) = {Pmf, AUm) ■ 

Notice that 

l(P./,A..)|<| l^ll^^^l < ^ .^.|^^^|. . 
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, > . Then we have 

We use part (iii) of Lemma (1) to obtain 

^^(«oll«mf + al\Aumf) + ^u{al\Aumf + al\A^/\mf) < 

< c||n„||(a^||„„|| + af |A3/2«„|)|An„|V2|A3/2^„|i/2 ^ 

< c\\um\\{alX^' + al)\A'/\mf/^\Aum\'/^ + K2 . 
By Young's inequahty we have 

^^("ollwrnf + al\Au^\^) + '^{al\Aura\^ + a?|A3/2u^|2) 

< c\\um\\\alX{^ + al)\ual)-Mum? + ^2 • (19) 
We integrate the above equation over (s,t) and use (17) and (18) to obtain: 

aolhm(i)f + al\Au^{t)\'^ < al\\u^{s)f + al\Au^{s)\'^ + 

2cfci2 

(^2)4^4 

Now, we integrate with respect to s over (0, t) and use (18) to get 
t{al\\um{t)\\'' + a\\Aum{t)\'') < ^{tKt+kx)+eK2 + 

+ (^|ff^("oAr^ + ain'^ + th] , (20) 

for all t > 0. 

For t > -v- we integrate with respect to s over the interval {t — ^i- , t) 



..^ .xx^^B^c.,.^ .v.^xx x^oi^^v^u o ^v^x ^xx^ xxx,.^xv<.x y. 



VA/ 2 j/Ai 



(21) 



From (20) and (21) we conclude: 

al\\um{t)f + al\Aum{t)\'' < k2{t) (22) 

for all t > 0, where ^(i) enjoys the following properties: 

(i) k2 (t) is finite for alH > . 

(ii) k2{t) is independent of m . 

(iii) If u'" e V, but m'" ^ ^(^)) then fc2(t) depends on i', f, ao and ai. Moreover, 
in this case limj__,o+ k2{t) = 00 . 

(iv) limsupj^oo k2{t) = < 00. 

Returning to (19) and integrating over the interval {t,t + T), for t > and r > 
and using (22) wo got 

ft+T 

{al\Aum{s)\^+al\A^/^Um{sW)ds < h{t,r) (23) 
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where ks{t,T) as a function of t satisfies properties (i)-(iii) as k2{t) above. Also, 
there exists Ti large enough, depends on (aQlu'^l^ + aiHu'"!!^), but independent of 
m, such that 

- / {al\Aum{s)f + al\A^/'^Um{s)f)ds < 2RI for all t > Ti . 
i Jo 

estimate (via the vorticity) 
Let us denote Vm = aoUm + aiAum and qm = V x Vm- The Galerkin system (15) 
is equivalent to 

+ vAVm - Pm{Um X Qm) = Pmf- 

Let us take Curl of the above equation, keeping in mind that we have periodic 
boundary conditions, to obtain 

^ + uAqm - V X [PmiUm X Qm)) = V X P^f ■ 

Notice that V • = and that PmQm = Qm- Let us take the inner product of the 
above equation with 

^^\Qmf + l^hmf - (V X {Pm{Um X qm)),qm) = (V X Pmf,qm) ■ 

We use the identity 

/ {V X (f)) ■ ^dx = I • (V X ^)dx, (24) 
Jq Jn 

to reach 

^^kmP + i^lkmf - {Pm{Um X g^), V X g„) = {Pmf,V X Qm) ■ 

Notice that Pm(V x g^) = V x therefore 

^^|9mP + i^|kmf = {Um X ^m, V X Qm) + (/, V X g„) , 

and upon applying (24) 

^^kmP + i^lkmlP = (V X (u„ X qm),qm) + (/, V X g„) . 

For every divergence-free function (j), and for every -0 wo have the identity 
V X (0 X V) = -((/> • + (V' • V)0 . 

As a result, we have 

^^kmP + i^lkmlP = -{{Um ' ^)q'm,qm) + {q-m ' ^Urn,qrn) + (/, V X g^) . 

Thanks to the identity (3) wc have {{um ■ ^)qm.,(Im) = 0. Now, we estimate the 
right hand side of the above to get: 

^^kmP + t^llgmf < cl|g™||L4||u™|| + \f\\\qm ■ 
We use the Sobolev inequality (6a) and Young's inequality to find 
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and we use Young's inequality again to obtain 



~\qrnf + ^\\qrnf<^\qrnf\\Urn\\' + l\ff. 



Let us denote Zm.{t) = u'^X\^^ + 19^(^)1^) then 



dt - "^^'y ^^x\'\ 

We use (17) to obtain 
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„_l/2)'^^ 

Zm{t) < Zm{s)e'' 



for every 0,s <t. From the definition of Zm we observe 

z^s) < cial\Au,nis)f + al\A^/^u„,is)f + v^X^^). 
Now, we integrate with respect to s over and use (23) to get 



^7 ^^ , 2\i/2 

^^3(2,2) + ^ V 



e 



2:m(t) < 

Here again k^i^ enjoys the properties (i)-(iii) of k^if), mentioned above. 



Remark 1. Notice that by establishing the estimate (25) for \qm\ one indeed is 
providing an upper bound for the _ff'^-norm of u„i. Similar estimates for the H^- 
norm of Um can be also obtained by considering first the Galerkin system (15) 

dv ~ 

+ vAVm + PmB{Um, Vm) = Pmf , 

taking the inner product with Avm , and then following a sequence of inequalities 
and estimates to achieve an upper bound for ||wm||- 

Let us now summarize our estimates. For any T > we have 

(i) From (17): 

II l|2 ^ ^1 II ||2 ^ I. 

|Pm||ioo([o,T];y) — 7^2 '^^ IP"»llL~([0,T];y) ^ "^1 

(ii) From (18) we have 

UT) „„ „„. „2 h{T) 



"mllL2([o,T];£>(A)) ^ °^ iPm lll,2([o,T],H) " 



< o- or \\v, 

(iii) From (22) 



2 ^(t) 2 

ll^»n|lz,<~([r,T];D(A)) ^ ^2 °^ II^TO|ll,«>([r,T];H) ^ ^2(t) 
"1 

for any r e (0, T], where k2{T) ^ oo as r — > 0"*". 
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Next, we establish uniform estimates, in m, for and ^g^. 
Recall (15) 

^'"m{t) = -PmB{Um, Vm) - ^AVm + Pmf ■ 

From the above estimates and part (v) of Lemma 1 we have 

IM l|2 ^ chiT) 

||^Wm||L2([o_T],D(A)') ^ ~ i 

and 

\\PmB{Um-,Vm)\\D(Ay < c\Um\^^'^\\Umf'^\Vm\ + km 

Consequently 



\v h( Ml 2 ^ ckik2{T) 

\^mi^\Um,Vmj\\L^m^T],D(Ay) ^ — 1/2 , 



Therefore 



~;lr\\\,'^{[0,T];D{A)') ^ ^(^) 



and in particular 



dt 

II dUm ||2 ^ fc(r) 

II^IIl=([0,T],H)<^, 

where fc(T) is a constant which depends on u, Ai, /, ao, a\ and T. 

By Aubin's Compactness Theorem (see, e.g., Constantin and Foias [1988] and 
Lions [1969]) we conclude that there is a subsequence Um'{t) such that 

Um' u{t) weakly in £^([0, T],D{A)) , 
Um' ^ u{t) strongly in L'^{[0,T],V), and 
Um'^u in C([0,T],if); 

or equivalently 

Vm' V weakly in L^qq, T], H) , 

Wm' ^ u strongly in L^([0,T],y), and 
in C([0,T],£)(A)') , 
where f is given in (9). 

Let us relabel Um' and Vm' by and Vm respectively. Let w G -D(A), then 
from (15) we have 

{Vm{t),w)+P {Vm{s),Aw)ds+ {B{Um{s),Vm{s), PmW)ds = 
J to J to 

(v„(io),iy) + {f,Pmw){t - to) , 

for all to,t € [0, T]. Since Vm ^ v weakly in i^([0, T]; _ff) then i'm(s) v{a) weakly 
in H, for every s £ [0,T]\E, where \E\ = 0. In particular, there is a subsequence 
of Vm, which we will also denote Vm, such that Vm{s) — > v{s) strongly in V and 
D{Ay for every s ^ E. 
Now, it is clear that 

rt rt 



lim / {vm{s),Aw)ds = / (t;(s), j4w)rfs 
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also that lim \PmAw — Aw\ = lim |m; — Wm\ = 0. On the other hand 

m — ^cxD m — >oo 

I [\BiuMvm{s)),Pmw) - {B{u{s), vis), w{s)) n^^y ds\ < 41) +7(^2) ^^(3) _ 

J to 

Im = I / {BiUmis),Vm{s)),PmW{s) - w{s)) D(^Ay ds\ 
Jto 

by part (v) of Lemma 1 we have 

4'^ < f\\\Um{s)\\\Vm{s)\ \PmAw - Aw\)ds , 

A-j^ 'J to 

applying Cauchy Schwarz inequality 

< ^ WUmisWds^ WVmisWds^ \Prr.Aw - Aw\ , 

and hence lim ijn^ — 0. 

m— J-oo 

Im =\ / {B{Um{s) -Um{s),Vm{s),w)D{Ayds\ . 
Jto 

Again thanks to part (v) of Lemma 1 

4^ < / \\Um{s) - U{S)\\ \Vm{s)\ \Aw\ds 

A]^ J to 

and by Cauchy-Schwarz 

4^ < ^ WUmis) - uisWds^ IVmis^ds^ \Aw\ . 

Since Vm bounded in L'^{[0,T];H) and Um ^ u in L^{[0,T],V) we conclude that 

Um = . 

Finally, 

ij^^ = \ {B{u,v-Vm),w)v'ds\ , 

Jto 

by virtue of part (v) in Lemma 1, and since Vm ^ v weakly in L'^{[0,T]; H), we 
obtain 

lim 4) = . 

m—^oo 

As a result of the above we have for every t^^t ^ [0, T]\E 

{v{t),w)+v I {v{s),Aw)ds -^r / {B{u{s),v{s),w)D{Ayds = 

J to Jto 

{v{tQ),w) + {f,w){t-to) , (26) 

for every w G D{A). Notice that since ||wm(i) ||l°°([o,t],v'') ^i, and since Vm{t) 
v{t) strongly in y for every f e [0,T]\i?, wehave ||i'(i)||LooQQ 7-]^y/) < ki. Moreover, 
because D{A) is dense in V, (26) implies that v{t) G C{[0,T];V') or equivalently 

«(t) eC([o,T],y). 

In particular, from (26) we conclude the existence of a regular solution for the 
system (8). 
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Uniqueness of regular solutions 

Next wc will show the continuous dependence of regular solutions on the initial 
data and, in particular, we show the uniqueness of regular solutions. 

Let u and u be any two solutions of equation (8) on the interval [0, T], with initial 
values u(0) = u™ and u{{)) = Ti™ respectively. Let us denote v = {a^u + alAu), 
V = {oqU + alAu), du = u — u, and hy Sv — v ~ v. Then from equation (8) we get: 

+ uAv + B{Su, v) + B{u, Sv) = . 

The above equation holds in ^^([o, T], D(^)'), since Su belongs to ^^([O, T], D{A)), 
the dual space of i^([0, T], D{Ay), wc use Lemma 1 to obtain 

{■^v, 5u)d(a)' + J^(q:oII^w|I^ + Oi\\A5u\^) + {B{u, 6v), 5u)d(a)' = . 

Notice that (■^, <5u)^(^), = i^(ag|(5u|2 + a?||5uf ), (see, e.g., Temam [1984], 
Chapter III, Lemma 1.2). As a result we have: 

~{al\5u\'' + a\\\5uf) + + al\A5u\'') + 

Now we use part (vi) of Lemma 1 to get 

\^{al\bu^ + a?||^wf ) + v{(^l\\^vf + a\\Abu^) < 
c{\ufl^\Au\^l''\bv\ \\bu\\\Au\ \bv\ 
and by Young's inequality we have: 

\jy^\H^ + "?P^if ) + ^{^IWH? + a\\A8u\^) < 
< ^{\\u\\ \Au\ \\Suf + \Au\^\Su\ \\6u\\) + 

+ ^{al\\6ur + al\ASu\') 
< ^\Au\\al\6u\'+al\\6uf) + 

+ '^{al\\5u\\'' + al\A5u\^). 

Hence, 

{al\5u{t)\^ + al\\5u{t)r) < (a^|<5«(0)|2 + a?||^n(0)f ) exp ( /* ^^^^rfs^ • 

Since u G -^^^([0, T],D{A)) we conclude the continuous dependence of the solutions 
of (8) on the initial data on any bounded interval [0, T]. In particular, we conclude 
the uniqueness of regular solutions. 

□ 

Remark 2. Following the techniques introduced in [19] (see also [17] and [26]) we 
can easily show that if the forcing term, /, in equations (8) belongs to a certain 
Gevrey class of regularity then the solutions of (8) will instantaneously belong to 
a similar Gevrey class of regularity. Specifically, in this situation the solution will 
become analytic in space and time. In particular, one can also provide uniform lower 
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bounds for the radii of analyticity (in space and in time) for the solutions that lie in 

the global attractor (sec section 4 for the existence of a compact finite dimensional 
global attractor.) As a result of this Gevrey regularity one can also show that the 
Galerkin approximating solutions, introduced earlier, converge exponentially fast 
in the wave number m, a,s m oo (see, e.g., [13], [20], and [25]). Furthermore, one 
can use this Gevrey result to establish rigorous estimates for the dissipative small 
scales in equation (8) (see, e.g., [14]). 



4. Estimating the dimension of the global attractor 

Let S{t) denote the semi-group of the solution operator to equation (8), i.e. 
u{t) = S{t)u™. It can be easily shown, from the proof of Theorem 3 and Rellich's 
Lemma (see [1]), that S{t) is a compact semi-group. Let us recall (see (13)) that the 
ball B\ = {u &V : \\u\\ <^}is an absorbing ball, in the space V. Consequently, 
the equation (8) has a nonempty compact global attractor 

A = C\s>o{^t>sS{t)Bi) 

(see, e.g., [11], [21] and [32]). 

In this section we employ the trace formula (see, e.g., [10], [11], and [32]) to 
estimate the Hausdorff and fractal (box counting) dimensions of the global attractor 
A in terms of the physical parameters of the equation (1). First, let us recall the 
Lieb-Thirring inequality 

Lemma 4 (The Lieb-Thirring inequality). Let {V'jIjLi an orthonormal set of 
functions in (H)'' = H (B ■ ■ ■ (B H. Then there exists a constant Clt, which depends 
k- times 

on k, but independent of N such that 

/ (^ ij, (x) ■ {x)f/^dx <CltY. (V^,- {x) : {x))dx . (27) 

Jn j^-^ j^-^ JQ. 



Next we will present a new technical Lemma which we will use in estimating the 
dimension of the global attractor. 



Lemma 5. Let {'■Pj}f^i CV be an orthonormal set with respect to the inner prod- 



net [•, •] which is defined in (2), i.e. 

[ipi, Vj] = al{^pi, ifj) + al{{(pi, (fj)) = 6i. 



Leti!j{x) = {aoipjix),aij^ipj{x),ai-^ipj{x),ai-£^ipj{x)y , and 

f'^i^) = J2j'=i{'Pj{x) ■ fji^))- Then, there exists a constant Cp, which is indepen- 
dent of N , such that 



. 1/2 



V{x)\\l^ < ^ iSlJ^i^M^) ■■ V^,(x))dx I . (28) 
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Proof. Let e R , j = 1, . . . , A'', to be chosen later. By Agmon's inequality (7) 

we have 

N N 

j=i i=i 

N N N N 

caol ^i'fi I II H ^i'^i II + II H ^i'^i H I ^J^'^i I ' 
j=i j=i j=i j=i 

and by Cauchy-Schwarz 

AT N 

j=l J=l 

1/2 / \ 1/2 

Af TV \ / A' N ^ 



<c I agiE^^-^^i'+"iiiE^^-^^ii' "oiiEs^^^ii'+"iiE^^-^^^- 

j=i i=i / V j=i j=i 



2 



< c 



AT JV 

j=i i=i 



1/2 / \ 1/2 / X 1/2 



h^Eii^if +"?Ei^^. 

u=i / V j=i j=i 



2 

'J I 



Since [i-pi, ipj] = Sij we have 

al\f2^j{A-'/\j)ix)\'+ai\f2CjVj{x)f 
j=i j=i 

1/2 

<^(E^I)hoEii^.f +"?Ei^^. 



. 1/2 



i=i \j=i 

Let i € {1, 2, 3} be fixed, we choose = ipji{x), from the above we have 

A A / A . 

«?(E4(^))' ^ c(E^?i(^)) E / ^ ^Mx))dx 

Now we sum over i, i = 1, 2, 3, to reach 
which concludes our proof. 



□ 



Theorem 6. The Hausdorjf and fractal dimensions of the global attractor of the 
viscous Camassa-Holm (NS—a) equations, dniA) and dp{A), respectively, satisfy: 

MA) < MA) < c^^{a"'(^f".Gn^^f"} , 
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where G = is the Grashoff number and, as before, ^ = min ■! , — ^ !• . 

Proof. The linearized (8) equation about a regular solution u{t) takes the form 

^6v + uA6v + B{6u,v)+B{u,Sv) =0 , (29) 

where v{t) = OQU + afAu and Sv = aQSu+alAdu. Notice that 6u evolves according 
to the equation 

^5u + uA5u + {ap + alA)~'^[B{5u, alu + alAu) + B{u, alSu + alASu) = , 

(30) 

which we write symbolically as 

^Su + T{t)Su = . 
at 

Let 6uj{0), for j = 1,... ,N, be a set of linearly independent vectors in V, 
and let Suj{t) be the corresponding solutions of (30) with initial value Suj{0), for 
j = 1,. . . ,N. We denote 

Tr,{t) = Trace (Pjv(i)r(t)|p^(,)^) , (31) 

where P/v(t)F = Il5vi{t) + 'RSv2{t) + • • • + R(5uAr(t), and PAr(t) is the orthogonal 
projector of V onto PN{t)V with respect to the inner product [•, •] given in (2). 

Let {'Pj{t)}jLi be an orthonormal basis, with respect to inner product [•,•], of 
the space PnV, i.e. [<^i, ipj] = Sij , i,j = l,... ,N. We set 

= (ao^,,a,^^,,a,^^,,a,^c,,)^ . 

Notice that {ipjjtpk) = Sjk , j,k = 1, . . . ,N. We set 

N 

i''^{x,t) = ^{i)j{x,t) ■ipj{x,t)) 

N N 

Notice that by the Lieb-Thirring inequality (27) 

/ {i>{x,t))^''/^dx<CLTQN{t) , 

where Qwit) := Ylf^i t) : Vil^j{x,t))dx. 

Let us denote 9j{x,t) = al(pj{x,t) + alA(pj{x,t), for j = 1, . . . ,N. Prom (31) 
we have 

N 

TN{t)^J2[Tit)^,{;t),^,{;t)] 

N N N 

= uY^[Aipj,ipj] +X(i3(</jj,i;),<pj) +X(i?(M,^j),(pj) , 
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and by virtue of (4) we have 

N N 



Observe that 



N 



N 



N 



N 



Thus 



= / (VV'j(a;,<) : yi^j{x,t))dx = QN{t) • 

TN{t) = uQN{t) + jN{t) , 



(32) 



where JN{t) '■= X^jLi(-S(u, ^j), yij). Let us now estimate JN{t)- Using (4) and (3) 
we have 

N 

JnH) = Y}^{u ■ V)ej,^j) + {{ipj ■ V)n, Bj)] 

N 

again by (3) 

JV 

JN{t) = ^[-al{{u-V)ipj,Aipj) + al{{ipj • V)u, (^j) + • V)u,A(^j) , 

integrating by parts and using (3) 

N 3 g g N 



j=l fc=l 

JV 3 



j=i fc=i j=i fc=i 



N 3 



Therefore, 



|Jjv(i)| < c / {Vu{x,t) : Vu{x,t)y/^^l;^{x,t)dx + 
Jn 



+ ca 



■/ 



N 



1/2 



dx 



where ip'^{x,t) = ^) ' -^s a result we have 



\JN{t)\ <c (Vu(x,t) : 'S7u{x,t)y/'^^l;^{x,t)dx + 
Jn 



fe,i=l 
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Thanks to (28) we have 

\JN{t)\ <c {Vu{x,t) : Vu{x,t)y^^ij'^{x,t)dx+ 
Jn 

/ 3 2 \ 

+C'J'QTit)J^ I J , (33) 

and by the Holder inequaUty we get 

|Jiv(i)l < c||V«|La/2 f / (V'(x,t))i°/3dx) ^ +cQ]i\t)\Au\( [ i;\x,t)dx^ ' 
\Jq J \Jq 

Since [(pi,ipj] = 6ij we have J^tl)'^{x,t)dx = N. Therefore, the above gives 

3 /5 

\JN{t)\ < c\\Vuhm ( [ mx,t)Y''/'dx) + cQ]l\t)\Au\N'/^ . 



Applying the Lieb-Thirring inequality (27) we obtain 

\JN{t)\ < C CLT||Vu||i5/2 (^J^i^^ji^^t) ■■ y^i{x,t))d: 

+ cQ^,i\t)\Au\N^I\ 

that is 

\JN{t)\ < c\\^u\\L./.Q'jl\t) + cQ'^\t)\Au\N^''' . 
Applying Young's inequality 



\JN{t)\ < -^\\S/u\\% + ^Q^(t) +c|^.x|4/3 
then using Holder's inequality 

\JN{t)\ < ^||V«l|^//||V^||i//+^Q^(t) + c|At.|4/3 
and by virtue of the Sobolev inequality (6c) we obtain 



2 \ 1/3 



^2x1/3 



iV2xV3 



\JN{t)\ < ^||w||V2||«f/4|^^|3/4 + |Q^(^) + ,|^„|4/3 

Using Young's inequality again we reach 
Substituting in (32) we get: 



1/2 



TN{t) >-QN{t) - 3/4 KIKlr + (^i\M') - 

t2\ 1/3 



-e,..r(^) 
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Now, we require N to be large enough such that 

1 

hminf- / Tjv(s)ds>0. (34) 

According to the trace formula (see [10], [11] or [32]) such an N will be an upper 
bound for the fractal and HausdorfF dimensions of the global attractor. Observe 
that from the asymptotic behavior of the eigenvalues of the operator A there is a 
constant cq such that 

Aj>coAi ,7-2/3 for j = l,2,... . 

Therefore, since QN{t) is the trace of the operator A restricted to some subspace 
of dimension N, we have 

JV 

Qwit) > 5^ Aj > cXiN^/^ . (35) 
Let us require iV to be large enough so that 

z/Ai7V5/3 > dim sup (^^ £ \Au(s)f/^ds^ (^^^ ' 

and 

v\rN^I^> I 3/4 l™^^P?^ f \\u{,s)fl\al\\u{s)f^<x\\Au{,s)\^)ds. 

V^l^a^ ai T^oo J Jo 

For such an N the inequality (34) is satisfied, and thus N provides an upper bound 
for the fractal and Hausdorfi^ dimensions of the global attractor. 
By Holder's inequality we have 

T / T \ 2/"^ 

limsup;^ / \Au{s)\^^^ds < limsup ( ^ [ \Au{s)f ] 
T->oo T Jq t^oo yT Jq J 

and thanks to (14) we get 

1 

limsup — / \Au(s)\^^^d^ ^ , „ 
T^oo T Jq y -faf 

Therefore, from the above, (13) and (14) wc have 

C / x2/3 / ^3/8- 

dH{A) < dF{A) < cmax I G*/^ (^:^^ j > G^/^ 



/^2„2),l/2\ 



which concludes our proof. 

□ 
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5. Numbers of degrees of freedom in turbulent flows 



An argument from the classical theory of turbulence [27] suggests that there are 
finitely many degrees of freedom in turbulent flows. Heuristic physical arguments 
are used to justify this assertion and to provide an estimate for this number of de- 
grees of freedom by dividing a typical length scale of the flow, L, by the Kolmogorov 
dissipation length scale and taking the third power in 3D. The resulting formula 
is usually expressed explicitly in terms of the mean rate of dissipation of energy 
and the kinematic viscosity. In analogy with this heuristic approach we will derive 
here an estimate for the "dissipation" length scale (i.e., what would correspond 
to the Kolmogorov length scale) for the viscous Camassa-Holm (NS— a) equations 
in terms of the mean rate of dissipation of "energy" and the kinematic viscosity. 
We will also show that the corresponding number of degrees of freedom is propor- 
tional to the dimension of the global attractor. This, in a sense, suggests that in 
the absence of boundary effects (e.g., in the case of periodic boundary conditions) 
the viscous Camassa Holm equations represent, very well, the averaged equation 
of motion of turbulent flows. Hence, one is tempted to use the viscous Camassa- 
Holm equations as a closure model for the Reynolds equations, which represent 
the ensemble-averaged Navicr Stokes equations. Indeed, this idea motivated our 
studies in [4], [5] and [6], and it also led to a physical derivation in [23] (see also [5]) 
of the viscous Camassa-Holm (NS— a) equations, in the inviscid case, as averaged 
equations. 

As before, we denote v 
take the form 



UqU + afAu, hence equation (8) and equation (10) 



dv - 

— + iyAv + B{u,v) = / 

u{0) = u*" 



(36) 



We define 



lim sup — 



{al\\u{s)f + al\Au{s)\'')ds 



(37) 



the mean rate of dissipation of "energy" , and 

e = sup e(w™) , 

the maximal mean rate of dissipation of energy on the attractor. From equa- 
tion (33), and since J iJj'^{x, t)dx = N, we have: 

\JN{t)\ <cJ^{Vu{x,t) ■.\7u{x,t))^/^iP^{x,t)dx + 

+cCl!^Q]^\t)\Au\N'/^ , 

and by Holder's inequality we have 

\JN{t)\ < cWWuheW^J^he/.+cQ'^^mulN'/^ 

Again by Holder's inequality and (6) we get 

\jN{t)\ <c\Au{t)\{j^i,\x,t)dxy^'^ {jjr{x,t)f/'dxf'+ 

+cQ]i\t)\Auit)\N'/^ . 
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Using the Lieb-Thirring inequality (27) and / 'tjj'^{x,t)dx = N we obtain 
\JN{t)\ < c\Au{t)\Ny'^Q](\t) + cQ](\t)\Au{t)\N'/\ 

and, hence 

\JN{t)\<c\Au{t)\N'/''Q]^\t) . 

After applying Young's inequality to the above and substituting in equation (32) 
we obtain 

Therefore, in order to satisfy (34), and based on the above, it suffices to choose N 
large enough so that for every trajectory u{t) on the global attractor A we have 

liminf Ci'^QNit) - ci/-i/3ivV9|AM(t)|4/3]rfi > . 

T->CX3 T Jq 2 

Therefore, such a large A'' is an upper bound for the dimension of the global attrac- 
tor. Based on (35) it suffices to require 

T->oo T Jq 

for every solution u'" e A, i.e., 

i/^/^AiA^s/^ > c sup (limsup^ / \Au{s)\'^''^ds\ . 
On the other hand, using Holder's inequality and (37) we have 



sup I lim 



a sup — / \Au{s)\^/^ds] < sup limsup(- / |AM(s)|2rfs 



2/3 



2/3 



< sup (limsup(— ^ / (aQ||u(s)||^ + ai|A'u(s)|^)rfs 

( ( \^^^ ( ~ \^^^ 

Therefore, every N large enough that 

3/4 



N>c[ (38) 

is an upper bound for the fractal dimension of the global attractor, and hence is an 
upper bound for the number of degrees of freedom in turbulent flows. 

We set the dissipation length scale, in analogy with the Kolmogorov dissipation 
length scale in the classical theory of turbulence, to be 

4 = 



Then equation (38) leades to the following 
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Theorem 7. The Hausdorff and fractal dimensions of the global attractor of the 
viscous Camassa-Holm (NS—a) equations, dniA) and dF{A), respectively, satisfy: 



This estimate for the number of degrees of freedom is consistent witli the conven- 
tional estimate a la Kolmogorov-Landau-Lifshitz [27]. In particular, the number 
of degrees of freedom scales as the cube of the ratio of the domain size divided by 
the Kolmogorov dissipation length scale (times a factor involving the fixed ai). 



We observed earlier that the system (8) reduces to the Navier-Stokes for ao = 1 
and ai = 0. In this section we will fix = 1 and investigate the convergence of 
the solutions of the system (8) as ai 0+, and relate the limit to the Navier- 
Stokes equations. We will be studying further the relation between solutions of the 
system (8) and the 3D Navier-Stokes equation in a subsequent work. 

Theorem 8. Let f G H, ?i™ G V and ao = 1. Let and Vo,^ ~ u^^ + afAua-^ 
denote the solution of the initial-value problem (8) (or equivalently (36)). Then 
there are subsequences u^j , u^j , and a function u such that as a\ — > 0+ we have: 
('^) '^a{ ~^ strongly in ^^^^([O, oo); H); 

(a) u, weakly in L^^^([0, cx)); V^); 

(Hi) for every T G (0, oo) and every w £ H we have {u^3^{t),w) — > {u{t),w) uni- 
formly on [0,T]; 

(iv) v^j — > u weakly in i|^^([0, oo); H) and strongly in i|^^([0, oo); V). 

Furthermore, u is a weak solution of the 3D Navier-Stokes equations with the ini- 
tial data u{0) = u™ (for the definition of weak solutions to the 3D Navier-Stokes 
equations see [11] and [31] 

Proof. Let T > be fixed. Prom the proof of Theorem 3 and by passing to the 
limit one can show that the estimates (17) and (18) also hold for the exact solution 
of the system (8). That is 



Jo 

This implies that there are subsequences {ui } and {vi }, and corresponding func- 
tions u and v such that: 




6. Convergence to the Navier-Stokes equations 



|u„,(i)|2+a2||u„,(f)f </ci 



and 





u 



weakly in 



L\[Q,T]-V), 



and 



{v^{}^v 



weakly in 



L\[Q,T]-H), 



as a[ — »■ 0+ 
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Next we will use the above estimates and equation (8) to show that 



dt 



dt 



f 

Jo 



dUai (t) 



dt 



D{Ay 



dt < K{T) 



(39) 



for some constant K which depends on T, but is independent of ai. Indeed, from 
equation (8) (or equivalently (36)) we have 

dUai 



dt 



+ uAua,+{I + ai^A)-'^B{uoc„Voc^) = {I + ai'^ f . 



Thus, 



_-^dUa^{t) 



dt 



< iy\u^,\ + \A-\l + ai^A)-^Biua„v^,)\ + \A-^f\ . 



In order to prove (39) we only need to find the proper estimate for 
\A-\l + a,^A)-'B{u^„v„,)\ < \A-'B{u^„v^,)\ . 
Applying part (v) of Lemma 1 we obtain 

\A~^B{Uai,Vai)\ < c(\UaA^^'^\\Uaj_\\^/'^\Va^\+ Xl^^'^\Va^\\\Ua^\\^ 



< 2cAi 



1/4, 



^CKl I II "CKl I 



< 2cA//^||uaJ|(|u„J+ai2|Au„J) . 

As a result of the above estimates we have 

\A-'B{u^„Vc,,)\^ < 8c2A-^/' ^ (^a,'\\u^^\\2){ai^\Auc,A^)) 

< 8c'A^'/'fci (||w„J|2 + ai2|^y^j2^ ^ 

and by integrating the above estimate over the interval [0, T] we have 

r \A-^Biua,{t),v„,{t))fdt< ^^Sc^A^'/^fci . 

Jo 

From all the above we conclude (39). 

By virtue of the above estimates and Aubin's compactness Theorem (see, e.g., 
[11], [28], or [31]) there exists a subsequence, which will also be labeled by {w^j}, 
that converges to u strongly in -^^([O, T]; H). Furthermore, since 

r \A-'^'Ki{t) - u^iim'dt = {air f \\n^{t)fdt < {a{f^ , 

we have that v^j u strongly in L^([0,T];y), as a\ 0+; and that v{t) = u{t) 
a.e. in [0,T]. 

As a result of these estimates one can extract subsequences, which will be also 
labeled by {m^j} and {f^j}, respectively, and show that as a\ 0^ 

B{u^{,v^^) ^ B{u,u) = B{u,u) weakly in L2([0, T]; i?(^)') , 

by following an approach similar to that used in the proof of Theorem 3. This 
finishes the proof of the Theorem. 

□ 
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